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ABSTRACT 

Context. Eight Accreting Millisecond X-ray Pulsars (AMXPs) are known to date. Although these systems are well studied at high 
energies, very little information is available for their optical/NIR counterparts. Up to now, only two of them, SAX J1808. 4-3658 and 
IGR J0029 1+5934, have a secure multi-band detection of their optical counterparts in quiescence. 

Aims. All these systems are transient Low-Mass X-ray Binaries. Optical and NIR observations carried out during quiescence give a 
unique opportunity to constrain the nature of the donor star and to investigate the origin of the observed quiescent luminosity at long 
wavelengths. In addition, optical observations can be fundamental as they ultimately allow us to estimate the compact object mass 
through mass function measurements. 

Methods. Using data obtained with the ESO- Very Large Telescope, we performed a deep optical and NIR photometric study of the 
fields of XTE 11814-338 and of the ultracompact systems XTE J0929-314 and XTE 11807-294 during quiescence in order to look 
for the presence of a variable counterpart. If suitable candidates were found, we also carried out optical spectroscopy. 
Results. We present here the first multi-band (VR) detection of the optical counterpart of XTE J1814— 338 in quiescence together with 
its optical spectrum. The optical light curve shows variability in both bands consistent with a sinusoidal modulation at the known 4.3 
hr orbital period and presents a puzzling decrease of the V-band flux around superior conjunction that may be interpreted as a partial 
eclipse. The marginal detection of the very faint counterpart of XTE J0929-3 14 and deep upper limits for the optical/NIR counterpart 
of XTE J1807-294 are also reported. We also briefly discuss the results reported in the literature for the optical/NIR counterpart of 
XTE J175 1-305. 

Conclusions. Our findings are consistent with AMXPs being systems containing an old, weakly magnetized neutron star, reactivated 
as a millisecond radio pulsar during quiescence which irradiates the low-mass companion star. The absence of type I X-ray bursts and 
of hydrogen and helium lines in outburst spectra of ultracompact (P orb < 1 hr) AMXPs suggests that the companion stars are likely 
evolved dwarf stars. 

Key words. 

1. Introduction ciate SXTs with LMXBs containin g an old, weakly magn etized, 

neutron star, (for a review see, e.g. lCampana et aL| [l998). In the 

Some Low Mass X-ray Binaries (LMXBs) exhibit sporadic out- following we will refer to these systems as Low-Mass X-ray 

burst activity while for most of their time they remain in a state Transients (LMXTs) 
of low-level activit y; these systems a re commonly referred to 

as X-ray transients Iwhite et al1ll984l) . Historically, the classifi- It has long been suspected that millisecond radio pulsars are 

cation of X-ray transients is made on the basis of their spectral the spun-up products of sustained mass transfer onto neutron 

hardness. The outbursts of Soft X-ray Transients (SXTs), charac- stars in LMXBs. According to this scenario, the companions 

terized by equivalent bremsstrahlung temperatures < 15 keV, are in the LMXBs evolve and transfer matter onto the neutron star 

often accompanied by a pronounced increase in the luminosity on a long time-scale, spinning it up to periods as short as a few 

of their (faint) opti cal counterparts and by the onset of type I X- ms. This model has gained strong support about ten years ago 

ray bu rst activity dMaraschi & Cavalierdl977tlWooslev & Taaml thanks to the launch of the Rossi X-ray Timing Explorer (RXTE) 

119761) . During quiescence, these transient systems are very faint satellite that discovered kilohertz quasi-periodic oscillations 

in X-rays (10 32 - 10 33 ergs" 1 ) and their optical luminosity drops (kHz QPOs) as well as coherent oscillatio ns during type-I 

by as much as 6-7 mag, giving a unique opportunity for the X-ray bursts in a number of LMXBs (e.g. IStrohmaver et al.1 

study of their companion stars. These properties clearly asso- E996). These periodicities are interpreted as the millisec- 
ond spin periods of weakly magnetic neutron stars. In April 

Send offprint requests to: P. D'Avanzo 1998 the first Accreting Millisecond X-ray Pulsar (AMXP) 

* Based on observations made with ESO Telescopes at the Paranal was discovered, with a 401 Hz (2.5 ms) X-ray pulsations 

Observatory under programmes ID 077.D-0677(D) and ID 079.D- (SAX J1808.4-3658; |Wijnands & van der Klis| |1998|). This 

0884(A) and observations made with the ESO telescopes obtained from conclusion was further strengthened by the discovery of addi- 

the ESO/ST-ECF Science Archive Facility. tional systems in the following years. Eight such systems are 



2 



P. D'Avanzo et al.: The optical counterparts of Accreting Millisecond X-Ray Pulsars during quiescence 



now k nown: SAX J1808.4-3 6 58 dW jjn ands & van der Klisl 
[19981: IChakrabartv & Morganl [l998b : XT E J 175 1-305 
dMarkwardt e t all |2002|); X TE J0929-314 dGallowav et al.l 
120021); XTE ~J1 807-294 (iMarkwardt'eFall l2003al): XTE 
338 dMarkwardt & Swankl 120031: IStrohmaver et al l 
IGR J0029 1+5934 dGallowav et all 120051): HETE 



J1814- 
120031) : 



J1900.1-2455 dMorgan et all 120051) and SWIFT J1756.9-2508 
dKrimm et al.ll2007H 1 L These findings directly confirmed evolu- 
tionary models that link the neutron stars of LMXBs to those of 
millisecond radio pulsars, the former being the progenitors of 
the latter. All these systems are LMXTs, have orbital periods 
in the range between 40 min and 4.3 hr and spin frequencies 
from 1.7 to 5.4 ms. These eight accreting millisecond pulsars 
are well studied at high energies, espe cially in X-rays , both 
during outburst and quiescence (see IWiinandsl 120061 for a 
review). On the other hand, with the significant exception of 
SAX J1808.4-3658 and IGR J00291+5934, their optical/NIR 
quiescent counterparts are only poorly known. The optical light 
curve of SAX J1808.4-3658 in outburst and quiescence shows 
variability modulated at the orbital perio d , in antiphase with 
the X-ray light cur ve dGiles et al.l 119991: iHomer et all |200U 
Campana et al. 2004]). This is unlike other quiescent transients 
that normally show a double-humped morphology, due to an 
ellipsoidal modulation, and indicates that the companion star is 
exposed to some irradiation. The same behaviour is seen in the 
quiescent optical light curve of IGR J00291+5934 (D'Avanzo 
et al. 2007fl In a theoretical study on SAX J1808.4-3658 
iBurderi et alj d2003l) proposed, on the base of previous works 
dStellaetal.ll 1994 Campana et al. 1998), that the irradiation 
is due to the release of rotational energy by the fast spinning 
neutron star, switched on, as a radio pulsar, during quiescence. 
Following this idea, Campana et al. (2004) and D'Avanzo et 
al. (2007) measured the required irradiating luminosity needed 
to match the optical flux of SAX J1808.4-3658 and IGR 
J0029 1+5934, respectively, and found that it is a factor of about 
100 larger than the observed quiescent X-ray luminosity for 
both systems. Neither accretion-driven X-rays nor the intrinsic 
luminosity of the secondary star or the disk can account for 
it. So, these authors conclude that the only source of energy 
available within these systems is the rotational energy of the 
neutron star, reactivated as a millisecond radio pulsar. 

Among the eight accretion-powered millisecond X-ray pul- 
sars known to date, four of them are in ultracompact systems, 
with orbital periods shorter than 60 minutes. Three of these ultra- 
compact systems, XTE J175 1-305, XTE J0929-314 and XTE 
J 1807 -294 are remarkably homogeneous, with measured orbital 
periods of 42.4, 43.6 and 40.1 minutes respectively, well below 
the minimum period for a system with a donor compose d pri- 
marily of hydrogen (f or b < 80 min; Rap paport et al.lll982l) . 

Optical and NIR observations performed in the past by dif- 
ferent groups only le d to deep upper lim its for the counterparts 
of XTE J 175 1-305 dJonker et all 120031) and XTE J1814-338 
dKrauss et al.ll2005|) or to th e detection of very faint candidates, 
if any (Monell ietalll2005l for XTE J0929-314). No observa- 
tions of XTE J 1807-294 during quiescence have been reported 
to date, while HETE J 1900. 1-2455 is still in outburst and SWIFT 



1 Very recently, episodes of coherent ms X-ray p ulsations were dis - 
covered in archival data of the LMXTs Aql X- l dCasella et afll2008h 
and SAX J1748.9-2021 jAltamirano et alj|2008l) . We will not consider 
these sources in the present paper. 

2 A recent paper of Jon ker et"aTl {2008) pointed out evidence for a 
change in the quiescent level of IGR J0029 1 +5934 and for strong flaring 
episodes. A slight phase difference in the folded light curve is also found 
with respect to that reported in D'Avanzo et al. (2007) 



J1756. 9-2508 has been discovered very recently (June 2007). 
The intrinsic faintness of the targets during quiescence, in com- 
bination with the large interstellar absorption and high stellar 
crowding of the relevant fields clearly explain the lack of detec- 
tions at optical wavelengths. 



2. Observations and data reduction 

Optical and NIR observations of the field of XTE J1814-338, 
XTE J0929-314 and XTE J 1807-294 were earned out with 
the ESO Very large Telescope (VLT), using FORS1, FORS2 
and ISAAC. All nights were clear, with seeing in the range 
0.5" - 1.0". For XTE J1814-338 we collected from the ESO 
archival VR images taken with the FORS2 camera on nights 
2004 May 20-21. A set of about 8 images of 5-6 min integra- 
tion were obtained every night which cover about 30 - 40% of 
the 4.3 hr orbital period at each filter. From the ESO archive we 
downloaded VRI images of the field of XTE J0929-314 taken 
with the FORS1 camera on nights 2003 December 19-20-21. 
The dataset consists of a set of 6 images of 5 min integration 
that cover more than twice per night the 43.6 min orbital period 
at each filter. We observed the field of XTE J1807-294 with 
our approved program ID 077.D-0677(D) in service mode with 
the ISAAC camera. Observations were carried out in /-band 
on 2006 May 5 and covered about 92% of the 40.1 min orbital 
period of the system. In addition, we collected from the ESO 
archive VRI images taken with the FORS2 camera on night 2004 
June 10 that cover about 30% -60% of the orbital period. We also 
obtained optical spectra of the counterpart candidates of XTE 
J1814-338 and XTE J1807-294 with the FORS1 camera on 
nights 2007 September 2-3 during our observing run approved 
under program ID 079.D-0884(A). The complete observing log 
is presented in Table Q] 

Image reduction was carried out following standard pro- 
cedures: subtraction of an averaged bias frame, division by 
a normalized flat frame; NIR frames were reduced using the 
ISAAC pipeline data reduction jitter, part of the ECLIPSE0 
package. Astrometry was performed using the USNOB1.C0 and 
the 2MASSncatalogues. PSF-photometry was made with the 
ESO-MIDA50 daophot task for all the objects in the field. The 
photometric calibration was done against Landolt standard stars 
for VRI filters and against the 2MASS catalog for NIR filters. In 
order to minimize any systematic effect, we performed differen- 
tial photometry with respect to a selection of local isolated and 
non-saturated reference stars. 

Our spectra of XTE J1814-338 and XTE J1807-294 were 
acquired with the grism 300V, covering the wavelength range 
4000-9000 A. We used a 1" slit, resulting in an effective res- 
olution of R - 440. The extraction of the spectrum was per- 
formed with the ESO-MIDAS software package. Wavelength 
and flux calibration of the spectra were achieved using helium- 
argon lamp and observing spectrophotometric stars. 

A short X-ray observation of XTE J18 14-338 was carried 
out with the XMM-Newton satellite. Details about the observa- 
tion and data analysis are given in Sec. 3.1. 



http ://archive. eso . org/ 

http://www.eso.org/projects/aot/eclipse/ 

http://www.nofs.navy.mil/data/fchpix/ 

http://www.ipac.caltech.edu/2mass/ 

http://www.eso.org/projects/esomidas/ 
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Source 




U 1 observation 


Exposure 


Seeing 


Instrument 


Filter/ 






(YYYmmdd) 


(s) 


(arcsec) 




Grism 


XTE J 18 14- 


-338 


20040521.27086 


10 x 60 s 


0.8" 


VLT/FORS2 


V 






20040521.27194 


9 x 60 s 


0.8" 


VLT/FORS2 


R 






20040522.26303 


10 x 60 s 


0.6" 


\ ti rr-i /I— /"\ I~l CO 

VLT/FORS2 


V 






20040522.25830 


9 x 60 s 


0.6" 


VLT/FORS2 


R 






20070904.07468 


4 x 1 800 s 


0.5" 


VLT/FORS 1 


300 V 


XTE J0929- 


-314 


20031220.27210 


6 x 300 s 


0.6" 


VLT/FORS1 


V 






20031220.27646 


6 x 300 s 


0.6" 


VLT/FORS 1 


R 






20031220.28083 


6 x 300 s 


0.6" 


VLT/FORS 1 


I 






20031221.29201 


6 x 300 s 


0.7" 


VLT/FORS 1 


V 






20031221.29636 


6 x 300 s 


0.7" 


VLT/FORS 1 


R 






O HA') 1 1 OAAH 1 

zUU31zzl.3UU/l 


6 x 300 s 


0.7 


VLl/rUKM 


I 






ZUUj izzz.zouzy 


A v "3AA c 
X JUU s 


u. / 


VT T/FHP Q 1 
V Ll/rUKo 1 


V 






20031222.26467 


6 x 300 s 


0.7" 


VLT/FORS 1 


R 






20031222.26904 


6 x 300 s 


0.7" 


VLT/FORS 1 


1 


XTEJ1807- 


-294 


20040610.16380 


9x60s 


0.7" 


VLT/FORS2 


V 






20040610.17390 


9x60s 


0.7" 


VLT/FORS2 


R 






20040610.18748 


16 x 60 s 


0.7" 


VLT/FORS2 


I 






20060505.19875 


30 x 3 x 20 s 


1.0" 


VLT/ISAAC 


J 






20070903.19999 


1 x 600 s 


0.9" 


VLT/FORS 1 


300 V 



3. XTE J1 81 4-338 

XTE J18 14-338 was discovered in outburst on 2003 June 
5 durin g monitoring of the centra l Galactic plane with the 
RXTE dMarkwardt & Swankl 120031) . The detection of coher- 
ent pulsations at frequency of 314.3 Hz (3.2 ms), modu- 
lated on an orbital period of 4.3 hr, makes this source the 
widest accreting mil l isecond X-ray pulsar known to date 
(IMarkwardt & Swankl 120031: IMarkwardt etail l2003bl) . Many 
X-ray bursts were observed from XTE J1814-338 with the 
RXTE satellite. As suming the Eddington li mit for the bright- 
est of those bursts Strohmayer et al. (2003) obtained a source 
distance of 8.0 +1.6 kpc. X-ray observations carried out with 
the Chandra satellite during outburst led to a precise source po- 
sition and showed a featureless X-ray spectrum which is best 
fit by an absorbed power law plus blackbody model, where the 
equivalent hyd rogen column densi ty (No) is comparable to the 
Galactic value (Kr auss et al.ll2.Q05h . An R ~ 18.3 candidate op- 
tical counterpart was identified within the small Chandra error 
box during the outburst and was no longer visible during quies- 
cence down to a limiting magnitude R > 23.3 (3<x c.l.; Krauss 
et al. 2005). Optical spectroscopy of the candidate carried out 
during the outburst revealed prominent hydrogen and helium 
emission lines, supporting its association with XTE J18 14-338 
dSteeghsll2003l) . An X-ray heated accretion disk model is able 
to account for the BVR magnitudes measured during the 2003 
outburst but this model is unable to account for the /-band data, 
which are systematically brighter than predicted (Krauss et al. 
2005). A similar "IR excess" with respect to the X-ray heated 
disk model, was observed in the accretion-powered millisec- 
ond pulsar SAX J1808.4-3658 (Wang et al. 2001; Greenhill et 
al. 2006) and suggests for both sources a significant contribu- 
tion from the companion star or from a synchrotron-emitting re- 
gion dWang et aljioOltlKrauss et alj2005l:lGreenhill et alj2006t 
iRussell et alJ 120071) . No radio observations of this source have 
been reporte d to da te. Using RXTE data taken during outburst, 
Papit to et alJ J2007) performed a precise timing analysis mea- 
suring a spin-down behaviour of the source (v 7 x 10" 14 Hz 

s _1 ) and used this value to estimate a magnetic field of 10 s - 10 9 
G. No quiescent X-ray observations of this system have been 
reported to date. 



3.1. XMM-Newton observations 

XMM-Newton observed XTE 1814-338 on Sep. 6 2005 for 41 
ks. The observation is plagued by strong proton flares and only 9 
ks (pn data) survive after filtering. A faint source can be detected 
at the known source position with a signal to noise ratio of 2.6. 
The source count rate is (4.9 + 1 .9) x 10~ 3 c s _I (1 <x confidence 
level). The source is not detected in the MOS1 and MOS2 ex- 
posures. Assuming a power law spectrum with T = 2 and the 
outburst column density (Krauss et al. 2005) we can estimate a 
0.5-10 keV unabsorbed flux of ~ 3 x 10~ 14 erg crrT 2 s , con- 
verting to a quiescent luminosity of ~ 2 x 10 32 erg s~' (assuming 
a distance from Earth of 8 kpc). Despite that this luminosity is 
somewhat higher than what usually observed in AMXPs, we re- 
mark that the uncertainties connected to this value are very large 
due to the paucity of counts as well as to the spectral parameters. 

3.2. Optical counterpart during quiescence 

In all our V and R frames we found a point-like object inside the 
Chandra error box at the following coordinates (J2000): R.A. 
= 18 A 13 m 39*.04, Dec = -33°46'22".3 (0.3" error). This posi- 
tion is coincident with the one of the optical counterpart of XTE 
J1814-338 detected by Krauss et al. (2005) during the 2003 out- 
burst. The object has V ~ 23.3 and R ~ 22.5. A finding chart is 
shown in Fig.Q] 

3.3. Optical light curve 

Once we identified the candidate, we searched for variability 
of this source in each filter. We performed PSF-photometry of 
our candidate and of a selected sample of bright, non-saturated, 
isolated stars assumed to be non-variable in all our VR frames. 
The result of such differential photometry is that our source is 
variable with a clear sinusoidal modulation at the 4.3 hr or- 
bital period with an average semiamplitude of about 0.4 mag 
(Fig. This unambiguously identifies the source as the opti- 
cal counterpart of XTE J18 14-338 and represents the first de- 
tection at optical wavelengths of this source during quiescence. 
In both V and 7?-band the light curve shows a single minimum 
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Fig. 1. Finding chart for XTE J1814-338 (30" x 30"). The de- 
tection of the source during outburst (left panel, from Rrauss et 
al. 2005) and /?-band VLT image taken during quiescence (right 
panel). The source is dimmer but still clearly visible. 

Table 2. Results of photometry of XTE J1814-338, all the val- 
ues are uncorrected for reddening. In column five are reported 
the reddening parameters used to correct our optical photome- 
try, computed assuming E(B — V) = 0.29 + 0.03 mag (see Sec. 
3.3 for details). 



Filter 




Mean ma; 


mitude 


Semi amplitude 


A, 




(A) 






(mag) 


(mag) 


V 


5270 


23.29 ± 


0.04 


0.52 ± 0.08 


0.87 + 0.03 


R 


6440 


22.52 ± 


0.03 


0.32 ± 0.08 


0.68 ± 0.03 



around phase 0, i.e. at superior conjunction (when the neutron 
star is behind the companion) and a maximum around phase 0.5 
(0.48 + 0.02 and 0.43 + 0.02 respectively, based on the precise 
X-ray ephemerides of Papitto et al. 2007), likely indicating that 
irradiation of the companion star plays a crucial role and it is 
reminiscent of the optical light curves of SAX J1808.4-3658 
and IGR J0029 1+5934 during quiescence. All the results of our 
phase-resolved multi-band photometry are reported in Table [2] 

It is interesting to note that in the V-band light curve, the 
points between phases 0.05 - 0.17 are well below the sinusoid 
best-fit of the data. We carefully checked the V-band images 
used to obtain the measurement, and it seems that they are not 
affected by problems (e.g., cosmic rays, spikes, bad pixels). In 
addition, the comparison stars are not affected by this decrease. 
However, there is no clear evidence of a similar effect in the 
corresponding points of the R -band light curve taken contem- 
poraneously (Fig. 0. A tentative explanation for this observed 
phenomenon will be presented and discussed in Sec. 3.5. 



22.5 



23 



23 



23.5 



24 - 



24.5 



0.2 


-0.2 
-0.4 
-0.6 



0.5 



1 

orbital phase 
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I t 

- A II 

1 1 1 1 


.i.ii 
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y 1 ' W : 
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Fig. 2. R and V-band light curves of XTE J1814-338. Phase 
zero is defined as superior conjunction (i.e. when the neutron 
star is behind the companion) and is computed on the X-ray 
ephemerides of Papitto et al. (2007). Two orbital phases are 
shown for clarity, the best sine-wave fit is also shown. The light 
curves for a field star of comparable brightness are plotted offset 
below (triangles). The residuals of the V-band light curve to the 
sinusoid model are shown in the lower panel. 



3.4. Optical spectrum 

Our spectra of XTE J18 14-338 were acquired on the night 2007 
Sept. 3, with the FORS1 camera of the ESO-VLT equipped with 
the grism 300V, covering the wavelength range 4000-9000 A. 
We calibrated them observing spectrophotometric standard stars 
and corrected the flux for interstellar absorption (see Sec. 3.5 
for details on reddening correction). We accounted for slit losses 
by matching our V and R photometry. With a simple rescaling 
(by a factor 1.2) we were able to correct the flux discrepancy 
between spectra and photometry, showing that the spectral flux 
calibration was robust. We then computed an average of all the 
available spectra, in order to increase the signal-to-noise ratio. In 
our averaged spectrum we can clearly detect Ha, H/3 and a pos- 
sible Hel /15875 emission lines, superposed on a blue continuum 
(Fig.©, with equivalent widths of (19.0 ± 1.1) A, (36.4 ± 3.6) 
A and (9.4 ± 1.2) A, respectively. These emission lines are also 



visible in each single spectrum. Unfortunately, the S/N and the 
spectral resolution are not high enough to determine if these lines 
have a double-horned profile (the typical signature of an accre- 
tion disk) or not. 

3.5. Origin of the quiescent optical emission 

Th e X-ray mass function ( fx(M) ~ 2.0 x 10~ 3 M o ) computed 
by lMarkwardt et alJ (l2003bl) . combined with a neutron star mass 
of 1.4 (2.0) M Q , implies for XTE J1814-338 a minimum com- 
panion mass of 0.17 (0.21) M p. Based on the harmonic p roper- 
ties of the burst oscillations, Bhattac haryya et all ((2005) were 
able to derive 90% confidence intervals of 26° < i < 50° for 
the inclination of XTE J18 14-338. A similar constraint was 
derived by Rrauss et al. (2005) from the absence of X-ray 
eclipses or dips and assuming an absolute magnitude of the sys- 
tem during quiescence of My < 13.2. With these constraints 
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Fig. 3. Average optical spectrum of XTE J 1 8 14-338 during qui- 
escence. We masked some bad-subtracted sky emission lines. 
The flux is normalized to continuum. Ha, HB and a possible Hel 
/15875 emission lines are visible. 



0.5 




0.5 1 1.5 

M (M SUN ) 



Fig. 4. Mass-color diagram for the optical counterpart of XTE 
Jl 814-338. Isochrones (from Baraffe et al. 1998) are plotted for 
low-mass stars with ages 0.1 (stars), 0.5 (squares), 5 (triangles) 
and 10 (diamonds) Gyr. Horizontal and vertical solid lines de- 
lineates our measure of the color and the derived estimate of 
the companion mass, respectively. Vertical dashed lines mark the 
mass estimates obtained by Krauss et al. (2005). 



on the system inclination the companion mass should be in the 
range (0.2 < Mi < 0.5) M Q . Information about the compan- 
ion mass can also be obtained from our color photometry. To 
compute the color of the source we first need to correct our 
magnitudes for interstellar absorption comp uted with the rela- 
tion N H /E(B - V) = 5.8 X 10 21 ctrr 2 mag" 1 dBohlin et al.ll 1978b 
and assuming N H = (1.67 + 0.17) x 10 21 cm" 2 (Krauss et al. 
2005). The resulting color excess is E(B - V) = 0.29 ± 0.03 
mag. Using a standard extinction curve from|Fitz£aMck|([l999) 
we obtained the reddening parameters for our filters (Table |2). 
We thus computed a mean unabsorbed (V - R) = 0.58 + 0.09 
color and compared it (Fig. |4]i with the theoretical mass-color 
diagr ams computed for solar metallicity low-mass stars reported 
in Baraffe ^t al.l (119981) . The result is that the secondary's mass 
should be 0.8 < M c < 1.0 M (lcr c.l.). Using the X-ray mass 
function, we can derive for these masses the inclination of the 
system 13° < i < 15°. Assuming a uniform a priori distribution 
of inclination angles, the inferred inclination would be low, and 
therefore improbable (P ~ 3%). In addition, these inclinations 
are at variance with the estimates of Bhattacharyya et al. (2005) 
and Krauss et al. (2005) described above. 

The discrepancies in the estimate of the companion's mass 
can find a natural explanation under the hypothesis that the sec- 
ondary star is subjected to strong irradiation from the compact 
object, as suggested by the optical light curve. Since during 
our observations XTE J18 14-338 was in quiescence, we try 
to investigate the possible causes o f its optical/NIR emis sion. 
Following Campana et al. (2004) and lD'Avanzo et al.l d2007l) we 
now attempt to account for the optical spectrum flux (corrected 
for interstellar absorption) with the simple model of an irradiated 
star with a black body spectrum (for the details of the modeling 
see IChakrabartvlll998l eqs. [8]-[13]). We fit the data by using 
the irradiating luminosity (L,„ ), the distance from Earth (D), the 
radius of the companion star (R c ) and the albedo of the star (77,) 
as free parameters. To fill the Roche lobe of a 4.3 hr binary sys- 
tem, a main-sequence companion star of XTE J18 14-338 should 
obey the mass- radius relation R c = O.28(M c /O.lM G ) I/3 fl (see, 
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Fig. 5. Optical spectrum (dots) and VR broadband photometry 
(triangles) of XTE J1814-338 during quiescence. Data are cor- 
rected for interstellar absorption. The plotted line represents the 
contribution of an irradiated star. See Sec. 3.5 for details. 

e.g.. lFrank et al.llT992h . For each possible radius R c we can thus 
derive the binary separation a = (G(M X + M c )(P 0[b ) 2 /(4n) 2 ) l/3 
which is also a parameter foreseen in the Chakrabarty (1998) 
model. We obtain an acceptable fit to all the data (see Fig.[5]and 
Tabled, with a reduced^ 2 = 0.9 (18 degrees of freedom, n.h.p. 
of 58%). 

The result of our fit is that the required irradiating luminosity 
is Li rr ~ 1 x 10 34 erg s _1 . Taking this value as a lower limit 
for the spin-down luminosity of a classical rotating magnetic 
dipole, we can estimate a neutron star's magnetic field greater 



6 P. D'Avanzo et al.: The optical counterparts of Accreting Millisecond X-Ray Pulsars during quiescence 

Table 3. Fit parameters used to model the XTE J1814— 338 optical spectrum (see Sec. 3.5 for details). Errors are at 90% confidence 
level. 

Model ~L^ r Distance ~R C rf, "x 1 

(ergs' 1 ) (kpc) (Ro) (cm) 

Star (1.02 ± 0.08) x 10 34 11.75 ±0.54 0.37 ± 0.02 0.16 + 0.07 - 09 

Disc + star (2.74 ± 0.03) x 10 34 10.25 ±0.10 0.34 ± 0.01 0.75 ±0.01 < 1.8 x 10 s < 0.8 




2000 5000 
wavelength (A) 



Fig. 6. Optical spectra (dots) and VR broadband photometry (tri- 
angles) of XTE J18 14-338 during quiescence. Data are cor- 
rected for interstellar absorption. The plotted line represents the 
contribution of an irradiated star plus a residual accretion disk. 
See Sec. 3.5 for details. 



than 9.7 X 10 7 Gauss. Such value is in agreement with the es- 
timate of Papitto et al. (2007). The required irradiating lumi- 
nosity is about two orders of magnitude larger than the qui- 
escent X-ray luminosity of XTE J1814-338 {L x ~ 10 32 erg 
s _1 ; see Sec. 3.1). Such discrepancy is reminiscent of that ob- 
served for SAX J1808.4-3658 (Bu rderi et al. 2003 ; Campana 
et al. 2004) and IGR J0029 1+5934 dD'Avanzo et alJl2007l) and 
can be explained with the presence in the system of a relativistic 
particle wind from an active pulsar which irradiates the com- 
panion star. As reported in Table [3] we were able to estimate 
R c = (0.37 + 0.02) R Q and consequently, using the mass-radius 
relation reported above, M c = (0.23 + 0.08) M . Such values are 
typical of a main-sequence M-type star, in agreement with the 
estimate made by Krauss et al. (2005) for this system. On the 
other hand, the (V - R) color we measured is consistent with 
a late G or K-type main-sequence star. In addition, the optical 
spectrum of XTE J18 14-338 is well fitted with a blackbody with 
J-max = (5106 + 24) A (Fig. [5), equivalent to a surface tempera- 
ture of T ~ 5675 K (typical of a G-K star). With a typical mass 
and radius comparable to the Sun, a G or K star cannot fit the 
Roche lobe of a binary system like XTE J1814-338. However, 
this discrepancy between the observed spectral type of the com- 
panion star and the geometrical limits of the system finds a natu- 
ral explanation if we assume that the relativistic particle wind of 
the radio pulsar, which irradiates the companion, would increase 
its surface temperature making it appear as an earlier spectral 
type star. 



The simple model we considered above does not predict the 
presence of a residual disk during quiescence. To check for a 
more realistic solution we tried to fit our data with the model of 
an irradiated star plus a disk by using the irradiating luminosity 
(Lj rr ), the distance from Earth (D), the radius (R c ) and the albedo 
of the companion star (77*) and the internal disk radius (/?,„) as 
free parameters. We fixed the X-ray albedo of the disk to 0.95 
(Chakrabarty 1998). For each possible radius /?,„ we assumed 
an outer disk radius of 0.3a (where a is the binary separation). 
We obtained an acceptable solution for values of the fitting pa- 
rameters similar to the model of the irradiated star alone (see 
Tableland Fig. [6j». We can provide just an upper limit for the 
inner disk radius, given that the reduced x 2 of the fit reaches a 
minimum for values of the inner disk radius that go beyond the 
level of the neutron star surface. Within the limits of our fit, the 
residual accretion disk contributes about 20% to the observed 
optical flu:xQ. Such scenario might shed light on the puzzling op- 
tical light curve we presented in Sec. 3.2 and in Fig. [2] A possi- 
ble way to explain this V flux decrease can be found assuming 
that the companion star or an intra-binary shock front (created 
by the interaction of the relativistic pulsar wind and matter out- 
flowing from the companion star and the disk) are eclipsing the 
disk. Since the disk is very blue (peaking at A < 1000 A), we 
would expect a marked decrease in the V band (~ 0.2 mag) and 
a minor decrease (~ 0.1) in the /?-bancQ, which might be hid- 
den in the error bars of our photometry. As reported in Sec. 3.2, 
the eclipse duration is of at least 10% of P or b (unfortunately we 
do not have data between orbital phase 0.80 - 1.05). We can 
also derive an upper limit of 0.3P or b on the eclipse duration as- 
suming that the eclipse profile is symmetric around phase 1.0 
(Fig. [2j. So, the size of the region producing the eclipse should 
be (2^a/P or b)0.1P ol -b < R e < (2ra/P or b)0.3P ol -b which means 
0.6a < R e < 1.9a, where a = (G(M X + M c )(P orb ) 2 /(4j:) z ) l/3 
is the binary separation. Using M x -1.4 M Q and M c = 0.2 M 
we obtain 1.0 R B ^ R e $ 3.0 R Q . This might indicate that the 
extended region is producing the eclipse, rather than a compact 
feature as also suggested by the broad decrease in the observed 
V-band flux. However, if the companion is eclipsing the disk, 
the system should be close to edge-on but, as noted by Krauss et 
al. (2005), no eclipses or dips are visible in the X-ray light curve 
(z < 77°), weakening the hypothesis of the companion eclipsing 
the disk. 

Alternatively, as shown in Campana et al. (2004) for SAX 
J1808. 4-3658, the intra-binary shock front (if present) is ex- 
pected to have a blue spectrum, partially contributing to the 
whole observed optical flux. If the companion star is eclipsing 
this shock front, then the eclipse should be visible also at low- 
inclinations, the shock front being spatially close to the compan- 
ion star. As in the previous scenario, such an eclipse should be 
more marked in the bluer bands, however only the lower limit 

8 We also tried to fit our data with a disk model, but we found an 
acceptable solution only for large values of the irradiation luminosity 
(~ 10 3S erg s~') and of the disk inner radius (~ 10 10 cm). 

9 Assuming that around phase zero the whole contibution of the disk 
(~ 20%) to the observed optical flux is eclipsed. 
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of our estimate of R e is compatible with the expected size of the 
companion star of XTE J1814-338. 

Finally, if the companion star is not fully convective (but 
see Podiadlovski 1991) and has a temperature gradient between 
the hemisphere facing the neutron star and the other one persis- 
tently hidden, this can explain the flux decrease without invoking 
eclipses. 

4. Ultracompact Accreting ms X-ray Pulsars 

Binary systems with P or b < 80 minutes can be formed 
in at least two possible channels as summarized, e.g., by 
iDelov e & Bildsten (2003). The first scenario starts from sta- 
ble mass transfer onto a neutron star from an e volved main- 
seque n ce star close to core hydro g en exhaustion dNelson et alj 
1986; Fedorova & Ergma 1989; Podsiadlowski et al. 2002; 
Nelson & Rappaport 200 3|), or a H e-burning star dSavoniie et alj 
1986; Tutukov &Fedorova| [T989). According to this scenario, 



Table 4. Results of PSF-photometry of XTE J0929-314, values 
in column three are uncorrected for reddening. In column four 
we list the reddening parameters used to correct our optical pho- 
tometry computed assuming E(B - V) — 0.15 + 0.02 mag. The 
absolute magnitudes are computed assuming a distance of 8 kpc 
(Wijnands et al. 2005; see Sec. 4.1.1 for details). 



the system goes through the 40 minute period range twice. 
Such a system will initially evolve to orbital periods comparable 
to the ultracompact AMXPs and can reach a period minimum 
of about 10 minutes. The resulting ultracompact binaries have 
donor masses of 0.1 - 0.2 M as they pass through P or b ~ 40 
minutes on their way toward the minimum period. At this stage, 
the donors have masses significantly greater with respect to those 
measured for the ultracompact AMXPs known to date. However, 
systems evolving through P or b ~ 40 minutes on the way out from 
their period minimum have masses of the order of 0.01 M , more 
in line with the measurements. By this time the donors have be- 
come partially degenerate low-mass, low-luminosity donor stars. 
On the contrary, the second channel involves evolution through 
the common-envelope phase during an unstable mass transfer 
from an intermediate mass (2-6 M ) to the neutron star. The 
donor and the neutron star spiral-in to shorter orbital periods un- 
til the envelope is expelled. The remaining core of the donor 
could be either an He or C/O white dwarf star with M ~ 0.0 1M 
dBildstenl 12001 IDelove & Bildsten! 120031) . The main difference 
between these two evolutionary models is that according to the 
first one a residual fraction of hydrogen still should be present 
in the companion star as the system reaches P or b ~ 40 minutes. 
Since the evolution of these systems is still unknown, a detec- 
tion of the companion star will allow us to differentiate between 
these different formation scenarios. 



4.1. XTEJ0929-314 

The X-ray transient XTE J0929-314 was discovered by the 
All-Sky Monitor (ASM) on the RXTE satellite in 2002 April 
(iRemillard & Smifhll2002l) . RXTE observations revealed persis- 
tent 185.1 Hz (5.4 ms) pulsations and a binary orbital period 
of 43.6 minutes. The measured X-ray mass function (fx(M) - 
2.7 x 1O~ 7 M ) gives a minimum companion mass of M c = 
0.008M G (i = 90°, M x = 1.4M ) and implies M c < 0.03M o 
(95% c.l.) for a uniform a priori distribution of inclination an- 
gles (Galloway et al. 2002). Chandra observations carried out 
during the 2002 outburst showed a featureless X-ray spectrum 
(0.5 - 8.3 keV) that is well fitted by a power-law + black- 
body model with Nh consist ent to Galaxy li ne-of-sight val- 
ues (0.8 - 1.0) x 10 21 cm" 2 dJuett et al.l l2003l). V ariable opti- 
cal dGreenhill et al.l l2002; Cacella 2002) and radio dRupen et alj 
2002) counterparts wer e detected during o utburst. Optical pho- 
tometry performed by iGiles et all d2005h during outburst re- 
vealed a sinusoidal modulation probably due to a combination 
of emission by a hotspot and an X-ray heated secondary star. 



Filter 


A c 


Mean magnitude 


A, 


Absolute mag 




(A) 




(mag) 




V 


5270 


28.2 ± 0.4 


0.46 ± 0.02 


13.2 ±0.4 


R 


6440 


27.1+0.3 


0.36 ± 0.02 


12.3 ±0.3 


I 


7980 


26.9 ± 0.4 


0.28 ± 0.02 


12.1 ±0.4 



These authors also report the evidence of an excess in the R and 
/ band in the BVRI broad-band spectrum (similar to what ob- 
served in SAX J1808.4-3658 and XTE J1814-338), that may 
be due to synchrotron emission in matter flowing out of the sys- 
tem. The optical spectrum of XTE J0929-314, carried out dur- 
ing outburst, shows a strong feature around 4640 A, probably 
due to C or N emission, and no clear evidence for H/He emis- 
sion or absorption lines (Nelem ans et alJ 12006). During quies- 
cence, XTE J0929-314 was barely detected by Chandra in the 
0.5-10 keV range. The quiescent X-ray spectrum could be fit- 
ted with a sim ple power-law model and no thermal component 
was detected dWiinands et alj|2005l) . Deep optical observations 
carried out during quiescence led to the detection of a very faint 
(R ~ 27.2) source at a position consistent with the X-ray and 
radio error circles (Monelli et al. 2005). 

4.1 .1 . Optical counterpart during quiescence 

We averaged all our V, R and / frames in order to obtain a deep 
image of the field of XTE J0929-314 for each filter. We found 
a faint, point-like object inside the precise radio error circle 
dRupen et alj|2002l) at the following coordinates (J2000): R.A. = 
09 ,, 29 m 20 i .19, Dec = -31°23'03".2 (0.3" error). This position 
is coincident with the one of the optical counterpart detected by 
Giles et al. (2005) during the 2002 outburst (Fig. [7]). The source 
is not visible in the single VLT frames, and no check for variabil- 
ity could be performed. The results of our PSF VRI photometry 
are shown in Table [4J we note that we only have a marginal de- 
tection of the source but our photometry is in agreement with the 
results reported by Monelli et al. (2005). 

We can use our detection to evaluate the possible nature of 
the optical counterpart of XTE J0929-314 in light of the pos- 
sible evolutionary states of ultracompact systems presented in 
Sec. 4. Given that the source is detected during quiescence, we 
will assume that, if present, the accretion disk contribution is 
low with respect to the total observed flux, as discussed for XTE 
J1814-338 (Sec. 3.5) and found for other accretion powered 
millisecond X-ray pulsars (Campana et al. 2004; D'Avanzo et 
al.2007). 

Using a column density 7.6 x 10 20 < N H < 1.0 X 10 21 cm -2 
(Juett et al. 2003) we estimate a color excess of E(B - V) = 
0. 15 +0.02 mag and, assuming for the system a distance of 8 kpc, 
we estimated the expected absolute magnitudes of the differ- 
ent possible Roche-lobe filling companion stars and compared 
them with our photometry of XTE J0929-314. Following simi- 
lar work by Jonker et al. (2003) on the AMXP XTE J 175 1-305, 
we considered, for the evolved main-sequence star scenario, the 
post-period-minimum model forXB 1832-330 (M = O.O26M , 
R = 0.055R G ), a binary system with orbital period close to 40 
min (Table 1 in Podsiadlowski, Rappaport & Pfahl 2002) and 
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Fig. 7. Finding chart of XTE J0929-314. The detection of the source during outburst (left panel, from Giles et al. 2005) and VLT 
image taken during quiescence (central panel). In the right panel we show a detail of the VLT image, with the radio position (Rupen, 
Dhawan & Mioduszewski 2002) marked by a cross. The optical counterpart of XTE J0929-3 14 (marked with solid lines) is dimmer 
but still visible. 
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Fig. 8. Absolute VRI magnitudes of the optical counterpart of 
XTE J0929-314 (dots) compared to the theoretical models for 
irradiated donor stars in ~ 40-min period ultracompact X-ray 
binaries (computed for different values of the irradiating lumi- 
nosity. See Sec. 4.1.1 for details). Solid lines: the post-period- 
minimum model for XB 1832-330 (Podsiadlowski, Rappaport 
& Pfahl 2002). Dashed lines: He-burning star (from Tutukov & 
Ferodova 1989). Dashed-dotted lines: the hot C/O white dwarf 
scenario (Deloye & Bildsten 2003). Dotted lines: the hot He 
white dwarf scenario (Bildsten 2002, Deloye & Bildsten 2003). 



the helium star donor model reported in Table III in Tutukov & 
Fedorova (1989; M = 0.015M o , R = 0.046# G ). For the white- 
dwarf scenario, we used the model of Deloye & Bildsten (2003) 
which estimates for XTE J0929-314 either a hot C/O star or a 
hot He star (M = 0.01 1M , R = O.O41fl and M = O.O19M , 
R = O.O5O/? , respectively; Figs. 5, 6 in Deloye & Bildsten 
2003). In light of the results obtained above for XTE J18 14-338 
and given the low number of data points, we considered the sim- 
ple model of an irradiated star with a blackbody spectrum (for 
details of the modeling see Chakrabarty 1998 and Sec. 3.5). As 
shown in Fig. [8] we see that for an irradiating luminosity equal 
to the X-ray quiescent one (Lx ~ 4.5 x 10 31 erg s _1 , assum- 



ing a distance of 8 kpc; Wijnands et al. 2005) none of the pro- 
posed models can account for the observed magnitudes of XTE 
J0929-3 14. To match our photometry, we need an irradiating lu- 
minosity higher than about a factor of 8 (L,„- ~ 3.5 x 10 32 erg 
s~'). Taking this value as a lower limit for the spin-down lumi- 
nosity of a classical rotating magnetic dipole, we can estimate a 
lower limit of 1.8 x 10 8 G for the neutron star's magnetic field 
(in agreement with the estimate of Monelli et al. 2005 and the 
upper limit of Wijnands et al. 2005). 



4.2. XTE J1 807-294 

XTE J 1807-294 was discovered on February 21, 2003 reveal- 
ing a frequency of 190.6 Hz (5.3 ms) by the RXTE satel- 
lite during a Galactic plane scan (Markwardt et al. 2003a). A 
follow-up Chandra observation allowed to obtain a precise po- 
sition (1" uncertainty, 90% c.l.) and to determine an orbital 
periodicity of 40.1 minutes, making it the shortest orbital pe- 
riod of all accreting millisecond X-ray pulsars known to date 
(Markwardt et al. 2003a). No type I X-ray bursts have been 
detected, so the dista nce of XTE J 1807-294 is still unknown. 
Campana et al.l (120031) observed the source during outburst with 
the XMM-Newton satellite at a luminosity level of about 2 x 10 36 
erg s _1 (0.5 - 10 keV, assuming a distance of 8 kpc) with an 
X-ray featureless spectrum well fitted by an absorbed (Nh ~ 
5 X 10 21 c m" 2 ) blackbody p lus hard Comptonization model like 
COMPTT (lTitarchuk|[T994l) . Kilohertz quasi-periodic oscillla- 
tions (kHz QPOs) have been detected in the X-ray flux of XTE 
J 1807-294 during outburst with a freque ncy separation con - 
sistent with the 191 Hz pulse frequency (Linares et al. 2005). 
From XMM-Newton outburst observations, Kirsch et al. (2004) 
improved the determination of binary orbital parameters and 
estimated an X-ray mass function of the system of fx(M) = 
1.6 x 1O~ 7 M , which implies (assuming a neutron star mass of 
1.4 M and i = 90°) a minimum companion mass of 0.007 M . 
A re-analysis of these data enabled Falanga et al. (2005) to con- 
strain the companion mass to M c < O.O22M , suggesting that 
is a very low mass (likely white) dwarf. INTEGRAL observa- 
tions of XTE J1807-294 in the 0.5-200 keV carried out during 
the 2003 outburst (simultaneously with the XMM-Newton and 
RXTE observations) revealed a spectrum consistent with a com- 
bination of thermal Comptonization and blackbody, as found at 
lower energy by Campana et al. (2003 ) but with a higher equiva- 
lent temperature (Falanga et al. 2005). During quiescence, XTE 
J 1807-294 was observed with XMM-Newton and the source was 
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Fig. 9. /-band (Le/f) and /-band finding chart for XTE 
J 1807-294. The candidate is indicated by solid bars. The black 
circle represents the Chandra error box. As reported in Sec. 
4.2.1, field star "A" is probably affecting our optical spectrum, 
due to its spatial closeness to the target. 



not detected down to a 3<x upper limit of L x < 4 x 10 31 erg 
s" 1 (0 .5-10 keV) assuming a distance of 8 kpc ( iCampana et al.l 
120051) . No detection of the counterpart of the system at other 
wavelengths has been reported so far. 

4.2.1 . A search for the optical counterpart during quiescence 

Our VRIJ imaging dataset of XTE J 1807-294 covers, respec- 
tively, about 32%, 32%, 57% and 90% of the 40.1 min orbital 
period of the system (see Table [TJ. An analysis of our VRI 
images revealed the presence of a source within the precise 
(1") Chandra error circle (Fig. |9]l at the following coordinates 
(J2000): R.A. = 18*06 m 59 s .77, Dec = -29°24'30".l (0.3" er- 
ror). As can be seen in Fig. [9] this source is barely visible in 
our averaged /-band image, obtained under worse seeing con- 
ditions, and no solution can be obtained with aperture or PSF- 
photometry in this band. VRI phase resolved PSF-photometry 
gave no evidence for variability. The source appears to be almost 
constant at a level of V ~ 22.1 and R ~ 21.4^ If any variability 
is present, it should be lower than 0.1 mag. Assuming a column 
density Nh ~ 5 x 10 21 (Campana et al. 2003) we estimated a 
color excess of E(B — V) = 0.79 + 0.05 mag and compute the 
unabsorbed (V - R) = 0.22 + 0.12 color, typical of an A5-F2 
main sequence star. 

To further investigate the nature of this source, we obtained 
on 2006 Sept. 3 an optical spectrum in the wavelength range 
4000-9000 A. We extracted the spectrum and found no evidence 
for Balmer or He emission lines (Fig. [TOl. The spectrum is al- 
most featureless, with the exception of Ha and HB absorption 
lines, superposed onto a blue continuum, suggesting an early- 
type star (in agreement with the results of our color photometry). 
The fact that the flux density of the spectrum is a factor of ~ 5 
brighter than our VR photometry (after correcting for interstellar 
reddening) shows that it is probably contaminated by a nearby 
bright star (star "A" in Fig. [9j». Among the possible scenarios of 
formation and evolution of ultracompact X-ray binaries that we 
have reported in Sec. 4.1, only an evolved main-sequence com- 
panion star is supposed to show hydrogen features in its optical 
spectrum but they should be detected together with prominent 



10 When calibrated with respect to Landolt standard stars, our /-band 
magnitudes of the candidate counterpart of XTE J 1807-294 appear to 
be suspiciously brighter than expected from the extrapolated (V - R) 
color. We found that this behaviour is common to other sources in the 
field and we conclude that our /-band photometry is very likely affected 
by some systematic error and, therefore, we report here only our VR 
results. 
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Fig. 10. Optical spectrum of candidate counterpart of XTE 
J1807-294 during quiescence (Fig. [9]). The position of typical 
optical emission lines of SXRTs are marked. 

Table 5. Upper limits (3<x) on the presence of a source at the X- 
ray position of XTE J 1807 -294. Values in column three are un- 
corrected for reddening. In column four the dereddening param- 
eters used to correct our optical photometry computed assuming 
E(B - V) — 0.79 + 0.05 mag are reported. The absolute magni- 
tudes are computed assuming a distance of 8 kpc (Campana et 
al. 2005; see Sec. 4.2. 1 for details). 



Filter 


X c 


Limiting magnitude 


A,, 


Absolute mag 




(A) 




(mag) 




V 


5270 


> 24.3 


2.38 ± 0.05 


> 7.4 


R 


6440 


> 24.2 


1.86 ± 0.05 


>7.8 


J 


12500 


> 19.7 


0.63 ± 0.05 


> 4.6 



helium lines (Nelemans, Jonker & Steeghs 2006). The detection 
of H absorption features together with the non-detection of He 
lines in the spectrum of our source and the absence of variability 
in the light curve suggests that our candidate is not the quiescent 
optical counterpart of XTE J 1807 -294 but more likely an early- 
type main-sequence star. We then used 3<x upper limits of our 
photometry (reported in Table |5]l to constrain the possible com- 
panion star of XTE J 1807-294. Assuming a distance of 8 kpc we 
converted our VRJ upper limits into absolute magnitudes (Table 

13. 

As done in Sec. 4. 1 . 1 for XTE J0929-3 14, we compared the 
results of our photometry to a set of spectral energy distribu- 
tions computed for the model of an irradiated companion star. 
All the models we considered for XTE J0929-314 (Fig.® can 
be considered valid for XTE J 1807 -294. In addition, we also 
considered the model of a cold C white dwarf (M = 0.01 1M , 
R = O.O41/? ; Deloye & Bildsten 2003; Falanga et al. 2005). 
We assumed an irradiating quiescent X-ray luminosity of about 
1 x 10 31 erg s" 1 (consistent with the non-detection of Campana 
et al. 2005) or a spin-down luminosity of 1 x 10 34 erg s _1 , as 
measured for the accretion-powered millisecond X-ray pulsars 
XTEJ1814-338(Sec.3.5), SAX J1808.4-3658 (Campana et al. 
2004) and IGR J00291+5934 (D'Avanzo et al. 2007). As shown 
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Fig. 11. VRJ 3cr upper limits for the optical counterpart of XTE 
J 1807 -294 (arrows) compared to the theoretical models for irra- 
diated donor stars in 40-min period ultracompact X-ray binaries 
(computed for different values of the irradiating luminosity. See 
Sec. 4.2.1 for details). Symbols used for the different models are 
the same as Fig. [8] with the exception of the model of a cold 
C white dwarf (Deloye & Bildsten 2003, Falanga et al. 2005), 
marked with a dashed-dot-dot-dot line. 



in Fig. QT| our upper limits do not constrain any of the models, 
even assuming a high irradiating luminosity. 

4.3. XTE J1 751 -305 

The transient accreting millisecond X-ray pulsar XTE 
J1751-305 was discovered on 2002 April during regular mon- 
itoring observations of the Galactic bulge region by the RXTE 
satellite. The pulsar has a spin frequency of 435.3 Hz (2.3 ms), 
a signature of orbital Doppler modulation on an orbital period 
of 42.4 minutes; no X-ray eclipses were detected (Markwardt 
et al. 2002). The precise determination of the orbital parameters 
enabled Markwardt et al. (2002) to derive the X-ray mass func- 
tion fx{M) ~ 1.3 X 10~ 6 M Q , which implies a minimum mass 
for the companion star of 0.013 - O.O17M depending on the 
mass of the neutron star. These authors derived a lower limit 
for the distance of the source of 7 kpc, by assuming a mass 
transfer rate from the companion equal to the one measured for 
SAX J1808.4-3658. XMM-Newton observations carried out dur- 
ing the 2002 outburst revealed an X-ray (0.5-10 keV) feature- 
less spectrum consistent with a simple absorbed blac kbody plus 
powe r - law model with Ng = 9.8 x 10 21 atoms cirT 2 dMiller et alj 
120031) . iGierlinski & Poutanenl d2005l) presented the results of an 
extensive study of all the X-ray data of XTE J1751-305 taken 
with the RXTE and XMM-Newton satellites during its 2002 out- 
burst. These authors show that the broad X-ray spectrum of 
XTE J1751-305 can be fitted with a model of two soft compo- 
nents (which originate from a cool accretion disk and a hot-spot 
on the neutron star surface) and a hard component due to ther- 
mal Comptonization in a shocked region in the accretion col- 
umn. During quiescence, XTE J1751-305 was observed with 
Chandra and was not detected down to a 0.5 - 10 keV luminos- 



Table 6. Upper limits (3<x) on the optical counterpart of XTE 
J 175 1-305; values in column three are uncorrected for redden- 
ing (from Jonker et al. 2003). In column four are reported the 
dereddening parameters we used to correct for interstellar ab- 
sorption computed assuming E(B - V) = 1.69 + 0.02 mag The 
absolute magnitudes are computed assuming a distance of 8 kpc 
(Wijnands et al. 2005; see Sec. 4.3.1 for details). 



Filter A c Limiting magnitude A,! Absolute mag 

(A) (mag) 

R 6500 > 23.1 4.03 ±0.02 > 4.6 

/ 8000 > 21.6 2.91 ±0.02 > 4.2 

Z 9625 > 20.6 2.12 ±0.02 > 4.0 

J 12500 > 19.6 1.35 ±0.02 > 3.7 



ity upper limit of (0.2 - 2) x 10 32 erg s assuming a distance of 
8 kpc (Wijnands et al. 2005). 

4.3.1. Constraints on the optical counterpart 

The X-ray observations of XTE J 175 1-305 place some strong 
constraints on the geometry of its companion star. The non- 
detection of X-ray eclipses or dips and the assumption that the 
mass transfer is driven by gravitational radiation (Markwardt 
et al. 2002) constrain the orbital inclination to be in the range 
30° - 85° and the companion mass to be 0.013 - O.O35M . 
Nelson & Rappaport (2003) computed a model for the evolu- 
tion of a main-sequence, hydrogen-rich companion star of XTE 
J 175 1-305 that can evolve to the observed short period of about 
42 minutes. According to this model, the current companion star 
should be a low-mass, low-luminosity non-degenerate star with 
a residual fraction of hydrogen. An alternative scenario is pre- 
sented by Deloye and Bildsten (2003) and foresees a hot white 
dwarf companion with an evolved (He, C/O) composition. The 
detection (e.g. via optical spectroscopy) of hydrogen would thus 
be a key discriminant for determining the validity of these mod- 
els. Unfortunately, there are no detections of the optical counter- 
part of XTE J 175 1-305 to date, nor in outburst or in quiescence. 
An observational campaign at optical and NIR bands placed up- 
per limits on the presence of a star at the precise Chandra posi- 
tion of XTE J175 1-305 (Jonker et al. 2003). 

We used 3cr upper limits of the Jonker et al. (2003) photom- 
etry (reported in Table |6]l to constrain the possible companion 
star of XTE J 1 75 1 -305 and converted them into absolute magni- 
tudes assuming a source distance of 8 kpc (Table |5). Following 
the procedure discussed in Sec. 4.1.1 for XTE J0929-314, we 
compared the results of Jonker et al. (2003) photometry to dif- 
ferent models of an irradiated companion star. All the models 
we have considered for XTE J0929-314 (Pig. [HJ can be consid- 
ered valid also for XTE J1751-305. We assumed an irradiating 
quiescent X-ray luminosity of about 1 x 10 erg s _1 (consistent 
with the non-detection of Wijnands et al. 2005) and a spin-down 
luminosity of 1 x 10 34 erg s erg s _1 . As shown in Fig. [12] 
these limits are not deep enough to constrain any of the pro- 
posed models for the companion star of this system, even for a 
high irradiating luminosity value. 

5. Discussion and Conclusions 

We have presented an extensive study of the optical counter- 
parts of Accreting Millisecond X-ray Pulsars during quiescence. 
The results of our observational campaigns (also presented in 
Campana et al. 2004 and D'Avanzo et al. 2007), together with 
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Fig. 12. RIZJ 3<x upper limits on the presence of a source at 
the X-ray position of XTE J 175 1-305 (arrows) compared to the 
theoretical models for irradiated donor stars in ~ 40-min period 
ultracompact X-ray binaries (computed for different values of 
the irradiating luminosity. See Sec. 4.3.1 for details). Symbols 
used for the different models are the same of Fig. [8] 



data collected from public archives, provided the first compre- 
hensive optical study of the AMXPs class during quiescence. As 
we discussed in Sec. 1, the eight accreting millisecond X-ray 
pulsars known to date can be divided into two subclasses using 
the orbital period as a discriminant parameter: compact (P or b > 1 
hr) and ultracompact (P or b < 1 hr) systems. The compactness of 
these systems implies very low-mass companion stars as inferred 
from the measured X-ray mass functions. So, it is reasonable to 
foresee a very low optical luminosity during quiescence, when 
the main contribution to the total optical flux is expected to come 
from the companion star. On the other hand, in such narrow sys- 
tems, irradiation from the compact object might play a key role 
in the observed luminosity, by heating the inner face of the com- 
panion star (this ene rgy is then redistributed through the star; 
Podsiadlowski T99lT) . We reported the first detection of the qui- 
escent optical/NIR counterparts of XTE J1814-338 (Sec. 3.2). 
Phase resolved photometry of this source revealed a highly sinu- 
soidal optical light curve, modulated at the orbital period with a 
minimum when the neutron star is seen behind the companion, 
with respect to the observer. This behaviour, reminiscent of what 
observed for SAX J1808.4-3658 (Homer et al. 2001; Campana 
et al. 2004) and IGR J00291+5934 (D'Avanzo et al. 2007), sug- 
gests a companion star heated by the compact object. This is at 
variance with the "classical" ellipsoidal modulation observed in 
wider LMXTs. The irradiating luminosity required to account 
for the observed optical flux is in all cases at least one order of 
magnitude greater than the observed quiescent X-ray luminos- 
ity. An explanation for this observed "optical excess" could be 
found if we assume that these systems host an active millisec- 
ond radio pulsar. According to the standard model of formation 
of millisecond pulsars these neutron stars have been recycled in 
interacting low-mass X-ray binaries, i.e. have been spun up by 
accretion after substantial decay of their surface magnetic fields. 
During quiescence (i.e. when the accretion rate onto the neutron 
star drops by orders of magnitude or even stops) the radio pulsar 



may turn on and start to lose rotational energy in the form of a 
relativistic particle wind that can irradiate the companion star. 
This wind eventually can be stopped by the pressure of the ma- 
terial flowing from the companion star generating a shock front. 
The presence of such an extended feature may eclipse the light 
coming from a residual accretion disk (if present during qui- 
escence), providing an explanation for the puzzling behaviour 
of the V-band light curve of XTE J1814— 338 we presented in 
Fig.|2]and discussed in Sec. 3.5. 

As for the compact AMXP systems, the detection of the ob- 
served optical flux during quiescence of the ultracompact system 
XTE J0929-314 can be accounted for only by invoking an irra- 
diating luminosity about one order of magnitude greater than the 
quiescent X-ray luminosity (Fig. 8). Taking our estimates of the 
irradiating luminosity as a lower limit for the spin-down lumi- 
nosity of a classical rotating magnetic dipole, we could estimate 
the neutron star's magnetic field for both systems. The derived 
limits are in agreement with the standard scenario of weakly 
magnetized (10 7 - 10 8 G) neutron stars to be the progenitors of 
millisecond radio pulsars (Sec. 3.5 and 4.1.1). We note that no 
radio pulsations were detected from XTE J0929-3 14 during qui- 
escence. Beaming factor and intrinsic low radio luminosity were 
suggested as the mo st likely explanations for this non-detection 
dlacolina et alj|2009h . XTE J 175 1-305 and XTE J 1807-294 lie 
both on the Galactic plane and are therefore heavily absorbed 
(we estimated Ay ~ 5.2 mag and Ay ~ 2.4 mag, respectively). 
Even assuming a high spin-down luminosity (L, rr = 1 x 10 34 erg 
s _1 ), the upper limits of the photometry are consistent with a non 
detection of the quiescent optical counterpart. 

As we discussed in Sec. 4, there are two main basic and 
competing scenarios for the formation of ultracompact binary 
systems. The first involves an evolved main-sequence or a He- 
burning companion star (both with a residual presence of hydro- 
gen) while the second predicts for the presence of a He or C/O 
white dwarf. We note that according to both scenarios the com- 
panion star of these ultracompact systems should be a very low- 
mass (M c ~ 0.01 M G ) low-luminosity star. As shown in Figures 
7, 10 and 11, with our data we cannot constrain any of the evolu- 
tionary models for ultracompact binary systems. However, as we 
reported in Sec. 4, the detection of hydrogen from one of these 
systems would be a key discriminant in determining the validity 
of the first scenario (which involves an evolved main-sequence 
or a He-burning star) with respect to the second one (white dwarf 
companion). Optical spectroscopy of XTE J0929-314 carried 
out during outburst did not reveal any clear evidence of hydro- 
gen or helium lines (Nelemans, Jonker and Steeghs 2006). In 
addition, we note that, at variance with what is observed for all 
the compact systems, none of the ultracompact AMXPs showed 
Type-I X-ray bursts. These events originate from explosive nu- 
clear ignition on the surface of neutron stars. At a critical point, 
degenerate hydrogen and/or helium burning ignites explosively, 
sudd enly heating up the entire surface, enough to emit stro ng X- 
rays ( iMaraschi & Cavalierdfl977t IWooslev & Taamlll976l) . The 
absence of Type-I bursts among ultracompact AMXPs finds a 
natural explanation if we assume that there is no hydrogen and/or 
He left in the donor star, as expected in the case of a C/O white 
dwarf companion star, suggesting that the other evolutionary 
channel for these systens is not populated. 

With quiescent unabsorbed magnitude in the range R ~ 
22-23, compact AMXPs are well within the spectroscopic capa- 
bilities of telescopes of the 8-meter class. An estimate of the ra- 
dial velocity of the companion star can lead to a determination of 
the optical mass function and, consequently, to place constraints 
on the mass of the two stars. Firm upper limits on the neutron 
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star mass could be fundamental to discriminate between different 
equations of state. Phase-resolved spectroscopy of these sources 
would be the natural continuation of the work presented in this 
paper. 

Acknowledgements. We thank the anonymous referee for her/his useful com- 
ments and suggestions. PDA and SC acknowledge the Italian Space Agency for 
financial support through the project ASI I/R/023/05. PDA wants to thank Prof. 
A. Treves for useful suggestion and discussion. PDA acknowledge H. Sana, S. 
Hubrig and L. Faunez for useful help given during observations at the ESO-VLT. 

References 

Altamirano, D., Casella, P., Patruno, A., Wijnands, R., & van der Klis, M. 2008, 
ApJ, 674, L45 

Baraffe, I., Chabrier, G., Allard, R, & Hauschildt, P. H. 1998, A&A, 337, 403 
Bhattacharyya, S., Strohmayer, T. E., Miller, M. C, & Markwardt, C. B. 2005, 

ApJ, 619, 483 
Bildsten, L. 2002, ApJ, 577, L27 

Bohlin, R. C, Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 

Burderi, L., Di Salvo, T., D'Antona, R, Robba, N. R., & Testa, V. 2003, A&A, 

404, L43 
Cacella, P. 2002, IAU Ore, 7893, 3 

Campana, S., Colpi, M., Mereghetti, S., Stella, L., & Tavani, M. 1998, 

A&A Rev., 8, 279 
Campana, S., D'Avanzo, P., Casares, J., et al. 2004, ApJ, 614, L49 
Campana, S., Ferrari, N., Stella, L., & Israel, G. L. 2005, A&A, 434, L9 
Campana, S., Ravasio, M., Israel, G. L., Mangano, V., & Belloni, T. 2003, ApJ, 

594, L39 

Casella, P., Altamirano, D., Patruno, A., Wijnands, R., & van der Klis, M. 2008, 

ApJ, 674, L41 
Chakrabarty, D. 1998, ApJ, 492, 342 
Chakrabarty, D. & Morgan, E. H. 1998, Nature, 394, 346 
D'Avanzo, P., Campana, S., Covino, S., et al. 2007, A&A, 472, 881 
Deloye, C. J. & Bildsten, L. 2003, ApJ, 598, 1217 

Falanga, M., Bonnet-Bidaud, J. M., Poutanen, J., et al. 2005, A&A, 436, 647 
Fedorova, A. V. & Ergma, E. V. 1989, Ap&SS, 151, 125 
Fitzpatrick, E. L. 1999, PASP, 111, 63 

Frank, J., King, A., & Raine, D. 1992, Accretion Power in Astrophysics 
(Accretion Power in Astrophysics, ISBN 0521408636, Cambridge University 
Press, 1992.) 

Galloway, D. K., Chakrabarty, D., Morgan, E. H., & Remillard, R. A. 2002, ApJ, 
576, L137 

Galloway, D. K, Markwardt, C. B., Morgan, E. H., Chakrabarty, D., & 

Strohmayer, T. E. 2005, ApJ, 622, L45 
Gierlihski, M. & Poutanen, J. 2005, MNRAS, 359, 1261 

Giles, A. B., Greenhill, J. G., Hill, K. M., & Sanders, E. 2005, MNRAS, 361, 
1180 

Giles, A. B., Hill, K. M., & Greenhill, J. G. 1999, MNRAS, 304, 47 
Greenhill, J. G., Giles, A. B., & Coutures, C. 2006, MNRAS, 370, 1303 
Greenhill, J. G., Giles, A. B., & Hill, K. M. 2002, IAU Circ, 7889, 1 
Homer, L., Charles, P. A., Chakrabarty, D., & van Zyl, L. 2001, MNRAS, 325, 
1471 

Iacolina, M. N, Burgay, M., Burderi, L., Possenti, A., & di Salvo, T. 2009, A&A, 
497, 445 

Jonker, P. G., Nelemans, G., Wang, Z., et al. 2003, MNRAS, 344, 201 
Jonker, P. G., Torres, M. A. P., & Steeghs, D. 2008, ApJ, 680, 615 
Juett, A. M., Galloway, D. K., & Chakrabarty, D. 2003, ApJ, 587, 754 
Kirsch, M. G. F, Mukerjee, K, Breitfellner, M. G, et al. 2004, A&A, 423, L9 
Krauss, M. I., Wang, Z., Dullighan, A., et al. 2005, ApJ, 627, 910 
Krimm, H. A., Markwardt, C. B., Deloye, C. J., et al. 2007, ApJ, 668, L147 
Linares, M., van der Klis, M., Altamirano, D., & Markwardt, C. B. 2005, ApJ, 
634, 1250 

Maraschi, L. & Cavaliere, A. 1977, in X-ray binaries and compact objects, p. 
127 - 128, 127-128 

Markwardt, C. B., Smith, E., & Swank, J. H. 2003a, The Astronomer's Telegram, 
122, 1 

Markwardt, C. B., Strohmayer, T. E., & Swank, J. H. 2003b, The Astronomer's 
Telegram, 164, 1 

Markwardt, C. B. & Swank, J. H. 2003, IAU Circ, 8144, 1 

Markwardt, C. B., Swank, J. H., Strohmayer, T. E., Zand, J. J. M. i., & Marshall, 
F. E. 2002, ApJ, 575, L21 

Miller, J. M., Wijnands, R., Mendez, M., et al. 2003, ApJ, 583, L99 

Monelli, M., Fiorentino, G., Burderi, L., et al. 2005, in American Institute 
of Physics Conference Series, Vol. 797, Interacting Binaries: Accretion, 
Evolution, and Outcomes, ed. L. Burderi, L. A. Antonelli, F. D'Antona, T. di 
Salvo, G. L. Israel, L. Piersanti, A. Tornambe, & O. Straniero, 565-568 



Morgan, E., Kaaret, P., & Vanderspek, R. 2005, The Astronomer's Telegram, 
523, 1 

Nelemans, G., Jonker, P. G., & Steeghs, D. 2006, MNRAS, 370, 255 
Nelson, L. A. & Rappaport, S. 2003, ApJ, 598, 431 
Nelson, L. A., Rappaport, S. A., & Joss, P. C. 1986, ApJ, 304, 231 
Papitto, A., di Salvo, T, Burderi, L., et al. 2007, MNRAS, 375, 971 
Podsiadlowski, P. 1991, Nature, 350, 136 

Podsiadlowski, P., Rappaport, S., & Pfahl, E. D. 2002, ApJ, 565, 1107 
Rappaport, S., Joss, P. C, & Webbink, R. F. 1982, ApJ, 254, 616 
Remillard, R. & Smith, D. 2002, The Astronomer's Telegram, 88, 1 
Rupen, M. P., Dhawan, V., & Mioduszewski, A. J. 2002, IAU Circ, 7893, 2 
Russell, D. M., Fender, R. P., & Jonker, P. G. 2007, MNRAS, 379, 1 108 
Savonije, G. J., de Kool, M., & van den Heuvel, E. P. J. 1986, A&A, 155, 51 
Steeghs, D. 2003, IAU Circ, 8155, 2 

Stella, L., Campana, S., Colpi, M., Mereghetti, S., & Tavani, M. 1994, ApJ, 423, 
L47 

Strohmayer, T. E., Markwardt, C. B., Swank, J. H., & in't Zand, J. 2003, ApJ, 
596, L67 

Strohmayer, T. E., Zhang, W., Swank, J. H., et al. 1996, ApJ, 469, L9+ 
Titarchuk, L. 1994, ApJ, 434, 570 

Tutukov, A. V. & Fedorova, A. V. 1989, Soviet Astronomy, 33, 606 

Wang, Z., Chakrabarty, D., Roche, P., et al. 2001, ApJ, 563, L61 

White, N. E., Kaluzienski, J. L., & Swank, J. H. 1984, in American Institute 

of Physics Conference Series, Vol. 115, American Institute of Physics 

Conference Series, ed. S. E. Woosley, 31 
Wijnands, R. 2006, in Trends in Pulsar Research, ed. J. A. Lowry, 53 
Wijnands, R., Homan, J., Heinke, C. O., Miller, J. M., & Lewin, W. H. G. 2005, 

ApJ, 619, 492 

Wijnands, R. & van der Klis, M. 1998, Nature, 394, 344 
Woosley, S. E. & Taam, R. E. 1976, Nature, 263, 101 



